Shear-Selected Clusters From the Deep Lens Survey III: Masses 

from Weak Lensing 



Alexandra Abate^, D. Wittman^, V. E. Margoniner^, S. L. Bridle^, Perry Gee^, J. Anthony 

Tyson^, Ian P. Dell'Antonio^ 

ABSTRACT 

We present weak lensing mass estimates of seven shear-selected galaxy cluster 
candidates from the Deep Lens Survey. The clusters were previously identified as 
mass peaks in convergence maps of 8.6 deg^ of R band imaging, and followed up 
with X-ray and spectroscopic confirmation, spanning a redshift range 0.19 — 0.68. 
Most clusters contained multiple X-ray peaks, yielding 17 total mass concentra- 
tions. In this paper, we constrain the masses of these X-ray sources with weak 
lensing, using photometric redshifts from the full set of BVRz' imaging to prop- 
erly weight background galaxies according to their lensing distance ratios. We fit 
both NFW and singular isothermal sphere profiles, and find that the results are 
insensitive to the assumed profile. We also show that the results do not depend 
significantly on the assumed prior on the position of the mass peak, but that this 
may become an issue in future larger samples. The inferred velocity dispersions 
for the extended X-ray sources range from 250 — 800 km s~^, with the exception 
of one source for which no lensing signal was found. This work further establishes 
shear selection as a viable technique for finding clusters, but also highlights some 
unresolved issues such as determination of the mass profile center without biasing 
the mass estimate, and fully accounting for line-of-sight projections. A follow-up 
paper will examine the mass-X-ray scaling relations of these clusters. 
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Introduction 



Clusters of galaxies are the largest virialized structures in the universe. Tracing the 
distribution of these structures as a function of cosmic time provides a measurement of the 



growth of structure an d is therefore a sensitive probe of dark energy (IHaiman et al.l 12001 



Battye &: Wellerl l2003l ) . This method of constraining dark energy requires accurate cluster 
mass estimates as well as unbiased samples of clusters selected over a broad range in mass and 
redshift (and with a well-defined selection function). There are tw o well-established selection 
techniques — optical d etection of the mem ber galaxies (IGal 1120081 ) and X-ray detection of the 
hot intr acluster gas (IRosati et al.l 120021) — and two new er techniques — weak gravitational 



lensin g (IWittman et al 



(SZE, ICarlstrom et al 



2001 



Hennawi fc Spergelll2005l ) and the Sunyaev-Zeld'dovich effect 



2002h . in which the cosmic microwave background is modified during 



its passage through the large potential wells associated with the intracluster medium. 

Each method has its strengths and weaknesses. Optical detection is capable of finding 
even poor clusters with relatively little observational resources, but its cluster selection 
function is the least well understood because the connection between mass (whose clustering 
is predictable) and galaxies is not well understood. X-ray selection is the least affected by 
line-of-sight projections because X-ray emissivity scales as the square of the local density, but 
by the same token dumpiness within the cluster changes the observable at a constant mass. 
SZE has the advantage of being nearly redshift-independent, but large SZE surveys are not 
yet available. Lensing selection is based most directly on mass, but it is observationally 
difficult because galaxies behind the cluster must be imaged well enough to infer the shear 
induced by the cluster. Whatever method is used to select the clusters, the next step is 
to infer cluster masses. Here the nonlensing methods depend on assumptions about the 
dynamical state of the cluster, and lensing may be useful for calibrating the masses of 
subsamples. However, this is observationally easier than determining the lensing masses of 
indiv idual clusters because subsamples may be stacked (jjohnston et al. 2007 ; Rozo et al. 
2008h . 



In this paper we determine lensing masses of shear-selected clusters, so hereafter we 
focus on the issues involved in shear selection and lensing mass determination. The weak 
lensing shear signal induced by a cluster is a function of its mass profile and its redshift 
relative to the background source galaxies being sheared, and is independent of the nature of 
the matter and of its dynamical or physical state. However weak lensing mass estimates do 
suffer from line-of-sight projections: structure along the line-of-sight could artificially boost 
the lensing signal or cause false positive detections. In addition poor knowledge of the source 
redshift distribution, as is often the case with only photometric redshift estimates, limits the 
accuracy of the mass estimate. We describe our choices in deahng with these issues as they 
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arise later in the paper. 

This paper is organized as follows. We describe the selection of and previous work on 
these clusters in §21 In §3], we summarize the data and its processing to the point of shear and 
photometric redshift catalogues. In ^ we describe the fitting procedure used to estimate 
the masses. We present the results in ^ and summarize and discuss the implications in 
^ Throughout this paper we assume a flat universe with = 0.3 and Hubble constant 
Ho = 70 km s"^ Mpc'^ 



2. Cluster Selection and Previous Work 



The Deep Lens Survey (DLS: IWittman et a. 



conducted with the Mosaic imagers (IMuller et al. 



2002) is a deep 20 deg^ BVRz' survey 



1998|) on the CTIO and KPNO 4-m tele- 



scopes, with nominal exposure times of 18 ks in R and 12 ks in the remaining bands. The 
clusters in this paper were shear-selected from the 8.6 deg^ which had at least 9 ks coverage 
in R band as of March 2002. Here we briefly summarize thei r selection and spectroscopic 
and X-ray foUowup. We refer readers to IWittman et al.l (120061 ) for much more detail. 



Wittman et al.l (120061 ) ranked the shear peaks and then folded multiple peaks within 



a 16' Chandra ACIS-I field of view into candidate fields for X-ray followup. The top eight 
candidate fields were observed for 20 ks each with ACIS-I in Chandra A0 4, with the exception 
of Abell 781, which had archival Chandra data. IWittman et al.l (120061 ) also obtained Keck 
and CTIO Hydra spectroscopy of likely cluster members, for redshift confirmation and to 
identify any line-of-sight projections. Redshifts were found to range from 0.19 to 0.68. Only 
one candidate (Candidate 5) was found to be a line-of-sight projection, on the basis of 
the spectroscopy and the absence of any detected X-ray emission. At the same time, most 
other candidates were found to contain multiple X-ray-emitting clumps, even where only one 
shear peak was present. In some cases, the additional extended X-ray sources correspond 
to extensions of the c onvergence map conto urs away from the main peak. See figures 8 
to 12 and 14 to 17 in IWittman et al.l (120061 ) where the multiband optical imaging of each 
candidate, with conve rgence maps and X-ray contours overlaid, were presented. The maps of 
Wittman et al.l (120061 ) were not well sui ted to resolving r nultip le peaks because of the heavy 
smoothing, as the kernel described by IWittman et al.l (120061 ) is approximately equal to a 
Gaussian with cr = 6.25' . In this paper, we use the shear data to fit a mass model to each 
extended X-ray source. This eliminates the need for smoothing, so we expect this process to 
yield mass estimates even where the original maps did not show a clear peak. We also use 
improved shear data as described below. 
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Two of the candidates have aheady been pubhshed in further detail 



Wittman et al. 




0.68. 


Sehffal et ah 



2003[) reported the spectroscopy and strong lensing arc of Candidate 8 at z 
( 120081 ) conducted a detailed X-ray and shear-fitting analysis of Candidate 1 at z = 0.30 (Abell 
781; note that the candidates are ranked by shear signal, so that Abell 781 was the top can- 
didate). Abell 781 has multiple clumps which were fit simultaneously so as not to bias the 
mass estimates for individual clumps. In this paper, we present simultaneous mass profile 
fits to all seventeen X-ray-emitting clumps in all the candidate fiel ds, as well as an over all fit 
to each candidate field. We repeat the fitting of Abell 781 because ISehgal et al.l (120081 ) used 
a completely different code and a slightly inferior version of the data, and we show that the 
two sets of results are consistent. 



3. Data 



The DLS includes five widely separated 2° by 2° fields, labeled Fl throug h F5. Fl and 



F2 ar e in the northern hemisphere and were observed with the Mosaic camera (IMuUer et al. 



19981 ) on the Kitt Peak Mayall 4-m telescope, while the remaining fields are in the south and 
were imaged with Mosaic II on the Cerro Tololo Blanco 4-m telescope. Each field is divided 
into a 3x3 grid of 40' by 40' subfields. The imaging data set for each subfield consists of 20 
dithered exposures in each of four bands [BVRz'), with each exposure lasting 900 sin R and 
600 s in the other bands. Observations were done in R when the seeing was be tter than 0.9" 
FWHM, and galaxy shapes were measured in R only. IWittman et al.l (120021 ) contains full 
details of the survey design. As mentioned above, the original candidate selection was done 
on the subfields (totaling 8.6 deg^) for which at least 9 ks (10 exposures) of R band imaging 
had been completed by March 2002. In this paper, we use the full imaging data set, so the 
data are deeper, and in some cases may extend to a larger distance from the candidate. 

For each filter and subfield, the 20 ex posures were co n abined into a single stacked im- 
age using custom software as described in IWittman et al.l (120061 ). The original images are 



8192x8192 pixels, with gaps between the eight CCDs and a pixel size which varies over the 
field but averages to about 0.26" . The stacked images are 10000 x 10000 pixels with a uniform 
pixel scale of 0.257" and with no gaps, cosmic rays, or satellite trails. The stacked point- 
spread function (PSF) size range s from 0.8-0.9" i n R a nd from 0.9-1.2" in the other filters, 
depending on field and subfield. IWittman et al.l (120061) contains rnany rn ore details on the 



stacki ng procedure. A minor improvement over the IWittman et al.l (120061 ) and ISehgal et al. 



( 120081 ) processing is the use of a Lanczos-windowed sine function for pixel interpolation. 
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rather than a truncated sine function]^ 



Cataloging and photom etric redshift estimation proceeded much as in ISehgal et al. 

( 120081 ) ■ We ran SExtractor (IBertin fc Arnoutslll996l ) on the stacked images in dual-image 



mode by detecting in R and measuri ng in each of BVRz'. We determined shapes of galaxies 
us ing the "VM" method desc ribed in lHeymans et al.l (120061 ). This is a partial implementation 
of iBernstein fc Jarvid (l2002l ) which does not give higher weight t o lower-ellipticity s o urces , 
so it should not lose accuracy in high-shear regions as does the full IBernstein fc Jarvid (120021 ) 



method. We apply a correction factor of to the VM shea r s, bas ed on performance in 



the blind analysis of simulations presented by iHeymans et al.l (12006! ) (corrected for stellar 
contamination which was present in those simulations but not in these data), and we fold a 
10% shear calibration systematic uncertainty into the cluster mass uncertainties. 



We used BPZ ( lBenftezll2000l ) t o estimate photometric redshifts , using the HDF prior and 
template tw eaking as de s cribed in iMargoniner fc WittmanI (12008! ) (based on the algorithm 
designed by !llbert et al.! (12006! )). When compared against 328 galaxies with spectroscopic 
redshifts in the NASA/IPAC Extragalactic Database {not used in the template tweaking), 
the rms photometric redshift error per galaxy is 0.047(1 + z), with a bias of 0.01 7(1 + z) and 
no catastrophic outliers over the redshift range 0.02 — 0.70 (ISehgal et al.l 12008! ). However, 
these are relatively bright galaxies, and we expect that performance is worse for the more 
typical faint galaxies used in the shear fitting. We address that concern in §4.21 below. 

For each subfield, we matched the shape and photometric redshift catalogs, eliminating 
galaxies which failed one or the other procedure, and then applied masks around bright 
stars. We also cut sources with extremely large or exactly zero errors on the measured 
ellipticity, indicating a problem with the ellipticity measurement; saturated sources; and 
sources with a photometric redshift greater than 1.6. We imposed the latter cut because 
photometric redshift performance is predicted to be poor after the 4000 Angstrom break is 
shifted through the z' filter, and we have no spectroscopy in that redshift range to evaluate 
how it degrades. Finally, we stitched together the nine subfields in each field, checking that 
object attributes in the overlap regions were consistent and that no objects were duplicated. 



^Lanczos interpolation is often used in popular packages such as SWarp, but the r ationale is not w idely 
known. Sine interpolation has the virtue of being exact for band- limited data (e .g Bracewelll 120001 ) . but 



it converges very slowly. Simple truncation, as was done for the IWittman et all (|2006l ) images, leads to 
artifacts; one can see this by considering the effects in Fourier space of a top-hat window in pixel space. 
The Lanczos-n window smoothly truncates the sine function at its nth zero by multiplying it by another 
sine function whose first zero occurs there. This also provides the useful property that the derivative of the 
interpolation kernel is zero at its edge. 
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4. Fitting procedure 

The observed reduced shear components, g" and are calculated from the measured 
ellipticities of the source galaxies ei and 62 by dividing by the shear responsivity correction 
Rg-. g° = Ci/Rg. We fit a model for the mass distribution of each cluster to the observed 
shear components, taking into account the full three-dimensional position of each source 
galaxy (r.a., dec, photo-z). 

We refer to each shear-selected cluster as a candidate and each X-ray detection within 
the vicinity of a candidate position as a "clump" . The X-ray positions give a better indication 
of a cluster center than the shear peak position because the X-ray centroids are more precise 
statistically. Chandra has 0.5" angular resolution, while the convergence maps have > 1' 
resolution. We number the candidates from 1 to 8 and label the X-ray clumps in their vicinity 
with lowercase letters. Where there are multiple shear peaks (which have X-ray confirmation) 
within what could be considered a single candidate, we label them with uppercase letters. 
For each clump the r.a. and dec. positions are given by the X-ray position, and the redshift 
is given by a spectroscopic measurement, except for clumps 2b, 2c and 7c. The photometric 
redshifts of these clumps are consistent with the spectroscopic redshifts of the main clump 
in their respective fields, so in these cases we use the main clump spectroscopic redshift. The 
information on the clus ters according to the l abelling described above is presented in Table 



dJ Note that Table 2 of IWittman et al.l (120061 ) contains a typographical error in the position 



of Candidate 1; the correct position for the candidate field is 09:20:50 +30:27:44 (J2000). 

In addition to the quality cuts mentioned in Section [3] we remove any sources which 
have redshifts less than Zc + 0.2, where is the redshift of the cluster. This is to avoid 
cluster member contamination in the sample of source galaxies. Cluster members scattering 
to higher redshifts would reduce the observed shear. When the fit applies to multiple clumps 
we take Zc to be the redshift of the highest redshift clump. The effect of the minimum 
photometric redshift cut on cluster member contamination is shown by Figure [1] for clump 
la; the same plot for the other clusters yields similar results. The number density as a 
function of radius is plotted before and after applying a minimum photometric redshift cut 
on the source catalogue. When there is no minimum photometric redshift (black stars) there 
is clear cluster member contamination at radii less than 10 arcminutes. When the minimum 
photometric redshift cut is applied (red circles) the number density is approximately flat 
with radius as shown by the fitted red dotted line. The gradient of the red dotted line is 
not significantly different from zero given the size of the scatter of the red circles. This 
implies there is no cluster member contamination, so no additional radius cut or correction 
is required. 



We use only sources that are within 30' of a profile center, partly for computational 
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Table 1: Top 8 shear-selected clusters from DLS. The clumps within each candidate are listed 

in order of decreasing X-ray flux. 



Cand. 


Shear peak info. 




Clump info. 








Field 


R.A.« 


Dec." 


Clump 


R.A.^ 


Dec.'' 


tF 




lA 


F2 


09:20:27 


+30:30:44 


a 


9:20:26.4 


+30:29:39 


0.302 


(1) 


IB 




09:20:50 


+30:27:44 


b 


9:20:53.0 


+30:28:00 


0.291 


(1) 


IC 




09:21:11 


+30:27:39 


c 


9:21:10.3 


+30:27:52 


0.427 


(1) 










d 


9:20:11.1 


+30:29:55 


0.302 


(1) 


2 


F3 


05:22:17 


-48:20:10 


a 


05:22:15.6 


-48 


18:17 


0.296 


(3) 










b 


05:21:59.6 


-48 


16:06 


0.296 


(2) 










c 


05:21:47.6 


-48 


21:24 


0.296 


(2) 










d 


05:22:46.6 


-48 


18:04 


0.210 


(5) 


3 


F4 


10:49:41 


-04:17:44 


a 


10:49:37.9 


-04 


17:29 


0.267 


(3) 










b 


10:49:50.7 


-04 


13:38 


0.068 


(4) 


4 


F4 


10:54:08 


-05:49:44 


a 


10:54:14.8 


-05 


48:50 


0.190 


(4) 


5 


F5 


14:02:12 


-10:28:14 










6 


F5 


14:02:03 


-10:19:44 


a 


14:01:59.7 


-10:23:02 


0.427 


(5) 


7 


F2 


09:16:00 


+29:31:34 


a 


09:15:51.8 


+29:36:37 


0.530 


(5) 










b 


09:16:01.1 


+29:27:50 


0.531 


(5) 










c 


09:15:54.4 


+29:33:16 


0.530 


(2) 


8 


F4 


10:55:12 


-05:03:43 


a 


10:55:10.1 


-05:04:14 


0.680 


(6) 










b 


10:55:35.6 


-04:59:31 


0.609 


(5) 



"Position of shear peak in DLS. 

''Position of X-ray peak in Chandra follow-up. 

'^Redshift sources: (1) Gcllcr ct al. 2005; (2) Wittman ct al. 2006, photometric; (3) Wittman et 
al. 2006: CTIO 4-m/Hydra; (4) CoUess et al. 2001: 2dF; (5) Wittman et al. 2006: Keck/LRIS; 
(6) Wittman et al. 2003 
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Fig. 1. — Cluster member contamination illustrated by the number density as a function of 
radius for clump la. The result is typical of the rest of the clusters. The number density is 
plotted before (black stars) and after (red circles) applying a minimum photometric redshift 
cut on the source catalogue. The minimum photometric redshift cut alone is adequate in 
removing cluster member contamination. 
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efficiency but mostly because sources at large radii will contain very little information on 



the c luster mass and be influenced by o ther nearby masses (see Fi gure 8 of I Johnston et al. 



20071 ) . Following iHoekstra et al.l (120021 ) and ICypriano et al.l (12004 ) we additionally remove 



sources within V of a cluster center. This is because in the most inner regions the mass 
densities could be closer to the critical value, leading to unmodeled strong lensing effects. 

For each model the best-fit parameters were obtained through evaluation of the 
statistic, defined as 

where g° and g2 are the measured shears. In Eq. [T]the i denotes the different source galaxies. 
The advantage of this approach is that it avoids data binning. We also use index m to 
distinguish between the two ellipticity components so is the predicted gi or g2 component 
for the source galaxy given a certain mass model. The error is the error on the ellipticity 
component which comes from the shape noise {(Tsn), the measurement noise (am) and the 
photometric redshift noise {a^z, see Section below) added in quadrature: = cr|^ + 
cr^ + 0"^^ . A separate measure of a^n for is used each galaxy. 

The shape noise is estimated by finding the standard deviation of gm values which 
have small measurement error (i.e. am < 0.002, leaving about 700 galaxies to compute the 
estimate). Then the shape noise asN is the average value estimated from gi and g2 and we 
find asN ~ 0.25. 

The predicted gl and g2 are calculated from the reduced shear: 

9 = T^ (2) 
1 — K 

where the form of 7 and k depends on the model (see Section U]T] below). Then the prediction 
for each shear component g{ and ^ff for a circularly symmetric distribution are: 

gP = -g cos 20 

92 = -fi'sin20 (3) 

where (f) is the angle subtended at the cluster center between the positive x-axis and the line 
from the cluster center to the position of the source galaxy. 



4.1. Mass Models 



We consider two different mass models, the singular isothermal sphere (SIS) and the 
Navarro, Frenk and White (NFW) profiles. The advantage of the SIS profile is that it has a 
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very simple form and only one free parameter, the velocity dispersion c^^,. The NFW profile 
has two free parameters, the concentration and M200, the mass within the radius at which 
the density is 200 times the critical density in the universe. 

The SIS profile is the density distribution found after solving the Boltzmann equations 
describing gravitational equilibrium of an isothermal gas. For the SIS profile the convergence 
and the shear are given by 

where 9e is the Einstein radius, which is related to the one-dimensional velocity dispersion 
of the isothermal sphere by 

Anal Dds . 
= — (5) 



s 



ay is the velocity dispersion (the mass proxy), c the speed of light, D^g is the angular diameter 
distance from the lensing mass to the source galaxies, and Dg is the angular diameter distance 
from the observer to the source galaxies. The shear from the SIS profile diverges at small 
projected radii, but the enclosed mass does not. With the 1' cut imposed for other reasons 
O, we do not attempt to fit these regions directly in any case and therefore do not encounter 
numerical stability problems. 

The NFW profil e is the distribution of dark matter predicted from N-body simulations 



(INavarro et al.l 119971 ). The characteristic radius of a cluster is defined as = r2oo/cc, and 
the projected radius relative to the center o f the lens as R. The anal ytic expression for the 
shear as a function oi x = R/vg is given by JWright fc Brainerdlboooh : 



^-^A{x) {x < 1) 

7(a:)= + 41n(i)] (a; = l) (6) 

r^l^B(x) (x > 1) 



wher e the functions A{x) and B{x) depend only on the dimensionless radius x (see lWright fc Brainerd 



2OOOL for their explicit form), 6c is a characteristic overdensity for the halo, pcr is the critical 
density of the universe at the redshift of the cluster and Scr is the critical surface mass 
density, 

D 

E,, = — (7) 

where is the angular diameter distance from the observer to the lensing mass. The 
convergence k is the ratio of the local value of the sur face mass density to the c ritical surface 



mass density, n = S(x)/Scr, where S(x) is given by (IWright fc Brainerd! l2000l ): 
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= 2rsScPcr 



VT- 



farctanh \/ 



(x < 1] 



ix = l) 



1 



, ^ , arctan , / f—!- 



fx > 11 



The dependence on the halo mass enters Eq. [S] and Eq. [HI mainly through the scale radius 
r^, and the dependence on the concentration mainly through and 5c- There is a weak 
dependance on cosmology through per, and x. 

We perform the fitting in two ways: 

(i) We fit one mass profile to each candidate. The high resolution of the X-ray data 
means that the X-ray position is a better estimator of the cluster center than the 
shear peak position. Therefore we vary the cluster center around the X-ray position; 
the cluster center is a free parameter. The X-ray position will usually be decentered 
slightly from the observed mass pe ak because of noise peaks superimposed on the mass 



distribution (see lClowe et al.ll2006l . and references within). Therefore even using a well 



resolved shear peak position could cause a systematic overestimate of the mass. If the 
offset is real, for example due to current or recent merging of clusters, then using the 
X-ray position would result in a systematic underestimate of the mass. At the same 
time, marginalizing over central positions over a large area biases the mass low. This is 
because fitting a profile positioned far from the true mass peak results in the prediction 
of a very low mass. We have confirmed this effect using simulated data. 

We therefore apply a prior on the cluster's X-ray/lensing offs et, and turn to the liter- 
ature for guidance on the size of the prior. ISmith et al.l (120051 ) examined this question 
observationally for ten massive clusters with high-resolution lensing maps from Hub- 
ble Space Telescope observations. In this case, the X-ray/lensing offset is much less 
affected by noise compared to the DLS clusters, which are less massive an d have lower- 
resolu tion lensing constraints due to lower background source densitv. ISmith et e 
( 2005 ) found offsets ranging from — 120 kpc, with an rms of 38 kpc. iKoester et al. 



(120071 ) examined 76 clusters identified in Sloan Digital Sky Survey photometric data 
and found a typical offset of 57 kpc between the X-ray centroid and the brightest 
cluster galaxy (BCG). For h = 0.7, the v alue assumed throu ghout this paper, this 
is roughly double the separation found by Smith et al. ( 20051). Assuming the B CG 
lies close to the bottom of the gravitational potential, the iKoester et al. (2007) re- 
sult should be more representative of our data because our clusters are of comparably 
modest richness. 
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We therefore adopt a prior of the form 

P = J\fe^r (9) 

where r is the angular separation on the sky between the highest flux X-ray peak 
position and the cluster center position, is the angular separation of a 57h^^ = 81 
kpc distance at the redshift of the cluster. We then marginalize over the cluster center 
positions. 

(ii) In the second approach, we simultaneously fit all the X-ray clumps in the vicinity of 
a candidate, to eliminate any influence one cluster may have on the inferred mass of 
another cluster. The shears from different clusters add linearly because they are small, 
so the predicted shear for the multi-profile fit is = + g^^^ + g^^ + summing 
up all the predicted reduced shears from each cluster, we use the exact X-ray peak 
positions as the profile centers in this fit. The g^^ etc are calculated as before (Eqs. 
Eland [3]) but the values of e.g. the distance ratios, and 16*1 (used to calculate g) will 
depend upon the clump in question. 



4.2. Photometric Redshift Errors 

Each of the source galaxies has a photometric redshift measurement and enters into the 
fitting separately, and each photometric redshift has an associated uncertainty that enters 
into the total uncertainty on the mass model. Here we calculate the relative magnitudes of the 
uncertainties from the photometric redshifts and from each galaxy's shear measurement. For 
this estimate of the relevance of photometric redshift errors only, we make the approximation 
that the observed (reduced) shear is approximately equal to the shear 7, 

^ = -^~7 (10) 

1 — K 

where k is the convergence. The shear can be written as 

7(^) = [ v{e-e')K{e')Se\ (ii) 



see 



Bridle et al.lll998l ). where V is the lensing kernel, and k is the convergence defined in 



Section 14. 1[ Using the definition of n and Eq. [7] one can write 



-Do 
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or 



'^GD,D,s ^^^^ 
where the mass model is now isolated on the left-hand side and the observables are isolated 
on the right-hand side. Uncertainties enter the right-hand side on an equal basis through the 
estimate of the shear 7 and through the distance ratio Dds/Ds which depends on the source 
photometric redshift (the distance to the lens is fixed through spectroscopic redshifts). In 
our data, the distance ratio uncertainty associated with a typical source galaxy's photometric 
redshift measurement is about 30% (expressed as a percentage of its distance ratio). The 
uncertainty in the shear estimate from any source galaxy is much larger, >500%, because a 
typical shear is 0.05 and shape noise — the uncertainty floor set by the random orientations 
of source galaxies — is 0.25. Therefore the per-galaxy distance ratio uncertainty of ~ 30% is 
negligible. 

Although we have shown that random photometric redshift errors are negligible on a 
per-galaxy basis, systematic errors in the photometric redshifts can still cause a bias on the 
measured cluster masses. To addr ess this question, we perf o rmed simulations of photometric 
redshift errors similar to those in iMargoniner fc WittmanI (120081 ). We ran multiple realiza- 
tions tailored to each cluster, starting with the actual DLS catalog for the subfield containing 
that cluster. For each realization, random photometric zeropoint errors were chosen accord- 
ing to the zeropoint uncertainties in each filter as derived from repeat observations. One of 
two priors for the distribution of re dshift and spec tral type as a func tion of magn i tude were 
chosen: either the HDF prior from iBenitezI (120001 ) . or the prior from lllbert et al.l (120061 ). In 
a given realization, for each galaxy, the R magnitude was used to choose a random type and 
redshift following the prior; these were used to generate synthetic colors and photometric 
noise and zeropoint errors were added. The photometric redshifts were then estimated using 



BPZ (IBenitez 



2000) , using either of the two priors at random, thus simulating the effect of 



uncertain knowledge of the prior. 



For a given lens, the induced shear scales linearly with the distance ratio as shown 
by Eq. [T^l Given the true and estimated redshift for each galaxy in a simulation, the ratio 
of true to estimated distance ratios was computed, and then averaged over the galaxies in 
the simulation. This yielded a distance ratio bias for that realization. The uncertainty in 
distance ratio bias was estimated by taking the rms of thirty realizations for each cluster. 
The typical bias was 5-10%, with the sign that the true distance ratio was smaller than 
estimated, which implies that the mass of the cluster is larger than estimated for a given 
shear. The typical uncertainty in the bias was 1-2%. 



We then use the distance ratio bias to "correct" the distance ratio calculated from the 
source photometric redshifts and the redshift of the cluster in question. The uncertainty 
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in this bias correction is negligible compared to that from other sources. A simple error 
propagation from the uncertainty in the correction to the mass estimator shows a ~1% 
effect on the mass estimate. The actual uncertainty must be larger becuase the simulation 
does not include some real effects, most notably a realistic variety of galaxy spectral energy 
distributions (SEDs). However, even if the uncertainty in the source redshift bias were 
doubled or tripled, the loosening of the mass constraints would be trivial, because other 
uncertainties dominate. At the same time, the 5-10% correction does capture the systematic 
effects of our filter gaps and photometry noise. 



5. Results 



For the NFW profile, we found that the lensing data provide very little constraint on the 
concentration parameter c. We therefore set the paramete r to c = 5, which is a typical value 



for the masses in question (ISeliakll200Cll: iDolag et al.ll2004r) . Setting the concentration to 6 (a 



high but reasonable value from iDolag et al.l |200J) instead of 5 changes the mass estimates 
by less than 10%, well within the uncertainty on the mass. When fitting the SIS profile, 
stability problems due to the divergence in the density at small radii were not an issue, due 
to the central 1' cut described above. To compare the results from the SIS fit and the NFW 
fit we convert the SIS into M200. They are related by 



200 



48 



2007rp, 



2G 



(14) 



where pcr is at the redshift of the cluster. 



We refer the reader to the figures presented in IWittman et al.l (120061 ) , which show the 
multiband optical imaging of each candidate, with the convergence maps and X-ray contours 
overlaid. 



5.1. Single Profile Fits 

The results from the single profile fits, described in Section WA\ item (i), are presented 
in Table [2] and Figured 

Table [2] shows the resulting masses with la errors after applying the prior for the cluster 
profile center position given in Eq. O We also include the offset in arcminutes between 
the best-fit centre position and the X-ray centre position; the average offset is about 0.4 
arcminutes. We find from simple simulations that a mis-center of this magnitude could still 
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bias the mass estimate by up to 5 per cent. Therefore not applying the prior (Eq. [9]) might 
have been significant for some of the clusters in our sample, given the size of their statistical 
errors. 

The mass found from the SIS profile is similar to the one found from the NFW profile; 
the two mass estimates agree to w ell within the la level . The mass found for candidate 8 is 
consistent with the estimate from IWittman et al.l (120031 ) once the shear calibration error in 



that paper has been corrected, see footnote 



mass estimates for Abell 781 (ISehgal et al 



3 in T able [2J Note that the previously published 



2008i) are not directly comparable to the new 



results presented in Table [2] because ISehgal et al.l (120081 ) performed simultaneous fitting, not 
sequential fitting of each clump as in Table [2l We therefore defer a comparison to §5.2[ in 
which we also fit all clumps simultaneously. 

Figure [2] presents the results from the single profile fits using the NFW profile. For 
each candidate M200 with its uncertainty is plotted at the cluster redshift. The uncertainty 
in the measured mass generally increases with redshift due to the reduction in the number 
of sources available for the analysis. Candidate 6 has a mass which is consistent with zero, 
and candidates 3 and 4 have masses which are nearly consistent with zero. This is a curious 
result for candidates which were selected on the basis of their shear. There are at least four 
effects at work here: 



• The uncertainties here include systematic and not just statistical effects. Thus, a 
candidate could be detected at, say 3a, but with systematic uncertainties its mass 
estimate may be closer to 2a above zero. 

• The fits here marginalize over a range of possible mass centers, whereas the origi- 
nal detections are precisely at the locations with the highest (smoothed) shear value. 
Noise fluctuations play some role in influencing this location, but at the same time, 
marginalizing over a large area will artificially lower the mass estimate and increase 
the uncertainty; see §4.11 for a discussion of these issues. Although we believe we have 
found a reasonable method for balancing these effects, the end result must be that 
the cluster mass estimate presented here has lower signal-to-noise than the original 
detection. 

• The original selection was made in the absence of photometric redshift information, 
while the mass estimates here use that information. Even in the absence of projections, 
this introduces competing effects: sources at ^phot > 1-6 were thrown away here, thus 
reducing signal-to- noise if those sources were truly at a reasonably high redshift. At the 
same time, signal-to-noise will be increased here to the extent that irrelevant sources 
(at redshifts less than the cluster redshift) can be discarded. 
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The use of photometric redshifts may correctly lower the signal-to-noise of projections 
which are widely enough separated in redshift. Even clusters which are real would have 
preferentially entered the sample if other groups of galaxies were projected near their 
line of sight. If the additional groups are at lower redshift, the cluster mass estimate 
will drop once source redshift information becomes available. If the additional groups 
are at higher redshift, it would be difficult to eliminate their effect with the simple 
redshift cut adopted here; they would have to be specifically identified and modeled 
out. The lowest-mass candidate in this paper. Candidate 6, in particular appears to 
have a foreground group of galaxies within a projected radius o f 2^ . Candidate 3 has a 
foreground group which is much more visible in the images of IWittman et al.l (120061 ) . 
and spectroscopically confirmed ai z = 0.069, but it is at a much larger projected 
distance (5'). A tomographic search making use of source redshift information pre- 
sumably would not have ranked these candidates so highly. We will present the results 
of a tomographic search in a future paper. 



5.2. Simultaneous Multiple-Profile Fits 

The results from the multi-profile fit, described in Section W7i\ item (ii), are shown in 
Table O Candidates 4 and 6 only had one X-ray peak so only one profile is fitted, but 
now the profile center position is set as exactly the X-ray peak position. The masses and 
la uncertainties in Table [3] are quoted after marginalization over the masses of any other 
clumps within the candidate field. 

Figure [3] shows the relationship between the masses from the single profile fit and multi- 
profile fit, with NFW results at left and SIS results at right. Comparing the two panels, 
we see that switching from single to multiple fits lowered the SIS mass estimates more 
than it lowered the NFW mass estimates. This suggests empirically that the NFW profile 
is less susceptible to shear from the neighboring clumps, a suggestion which is supported 
mathematically by the greater steepness of the NFW profile at larger radii. 

We now examine the effect of single vs. multiple fits on individual clusters. Candidates 
1 and 2, with the most X-ray-emitting clumps (four), are the two with significantly lower 
masses in the multi-profile fit as compared to the single-profile fit. This is simply due to the 
shear of the neighboring clumps contributing "mass" to the single-profile fit. Candidates 7 
and 8, with three and two X-ray-emitting clumps respectively, also have lower masses with 
the simultaneous-fitting approach, but this effect is not significant given their larger uncer- 
tainties. For the remainder, the mass estimates did not decrease with the simultaneous-fitting 
approach; of these, only Candidate 3 had more than one X-ray-emitting clump. Collectively, 
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Table 2: Results from single 



profile fitting 



Cluster 


z 


Centre offset 


NFW Result 


SIS result 


Published result 


No. 






(arcminutcs) 


M200 {10"MQ/h) 


a-jj (km/s) 


A/200 (lO"M0/h) 


M20O (lO^^Mg/h) 


sources 


lA 


0.302 


0.52 


4 2+0-6 


900t^^ 


4 1+0-5 


44+1.6a 


49411 


IB 


0.291 


0.34 


4 fi+°-^ 

^■"-0.5 


925+g 


4 t:+0-6 
+ •"3-0.5 


Q 7+2. la 
-J- ' -1.5 


50456 


IC 


0.427 


0.40 


c; 4+0.8 


1025t^9 


c fi+0-8 


/I 0+2. 7a 


37648 


2 


0.296 


0.31 


1 0+0.4 


700t^4 


1 q+0-4 




45249 


3 


0.267 


0.43 


U.D_Q 2 


400tl«9 


^•+-0.2 




42743 


4 


0.190 


0.48 


n 4+°-2 
^•^-0.2 


99r+85 


^■^-0.1 




40991 


6 


0.427 


0.49 


n 4+°-2 
^■^-0.4 


300t^J^ 


"•^-0.1 




28626 


7 


0.530 


0.31 


9 9+1.0 
^■^-1.2 


800t-° 


9 ^+1.3 




16354 


8 


0.680 


0.25 




79r + 121 


1 7+0.9 

-■^■'-0.7 


2.8tl.O'' 


19207 



"Converted from M500 v alue in ISehgal et al.l ll2008h . 

''Converted from a-u in IWittman et al .1 l|2003l '). Note: We found there was an error in the shear cahbration of 
IWittman et al.l ll2003h . With proper calibration, the weak lensing mass estimates and statistical uncertainties in that 
paper should be divided by a factor of 1.6. The mass presented here has had this correction applied. 




Fig. 2. — Mass from the single NFW profile fit as a function of the cluster's redshift. See 
the legend to relate the data points to the candidates. 
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1 2 3_, 4 5 6 1 2 3 _, 4 5 6 

single profile fit Mjuu (h solar masses) x1o" single profile fit M^uu (h solar masses) x1o" 



Fig. 3. — Comparison of masses from the single profile fit and the multi-profile fit. The 
left hand panel shows the results from the NFW profile and the right hand panel shows the 
results from the SIS profile. See the legend to relate the data points to the candidates. Note 
that for clarity the candidate 2 result from the NFW profile fit (red/dark star point, left 
hand panel) has been offset to the left by 0.1x10^^ Mq/H, which is much less than the size of 
its errorbar. 




Fig. 4. — Comparison of masses found from fitting the SIS profile and the NFW profile. The 
green (light) triangle points are the results from the multi-profile fit. The blue (dark) circle 
points are the results from the single profile fit. 
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Table 3: Results from multi-profile fitting 



Clump 


Peak 


z 


NFW 




SIS 


Published 


No. 




Value 




A/200 (lO^Mg/h) 


(Ti, (km/s) 


A-f200 (10"Mo/h) 


Af200{10"MQ//i) 


sources 


la 
lb 
Ic 
Id 


0.0188 


0.302 
0.291 
0.427 
0.302 


q c+0.4 

q Q + 0.6 

O.U_g 4 

2 0+° '^ 

0+"-^ 


750+3? 
625+^7 
725i?^ 
300t«°„ 


2 4+0.3 

2 0+°-6 
2+°-i 


44+1.6a 
4 q+2.7a 


48228 


2a 
2b 
2c 
2d 


0.0151 


0.296 
0.296 
0.296 
0.210 


-•^■"-0.4 

-■■■^-0.4 

1 c;+0.5 
-'^■^-0.6 


400+^^ 
650t?« 

^'3-126 


Q 7+0.3 
'-'•'-0.2 

^•"-0.5 

i+°-2 




51738 


3a 
3b 


0.0136 


0.267 
0.068 


U.O_g 2 



^■^-0.2 


400+p 
250 127 


4+"-^ 




44680 


4a 




0.190 


4+°-2 


375169 


0.31^2 




41085 


6a 




0.426 




^■^-0.3 


400i-° 


3+"-'' 




28632 


7a 
7b 
7c 


0.0119 


0.530 
0.531 
0.530 


1 6+1" 

1+°-^ 

^•■■--o.o 

^■^-1.0 


675+-^ 
n+215 

77^+117 

' '3-159 


1 5+^-" 

2.3tl:l 




19357 


8a 
8b 


0.0119 


0.680 
0.609 


1 3+"'^ 

2 4+0.9 


"J' 3-167 

800i?^, 


0.9t°j 

2.4t°:^ 


2.8il.O'' 


20289 



"^Converted from M500 value in ISehgal et al.) tOO^ . 
Converted from (T„ in IWittman et al.l l|2003i v see footnote b in Tabled 
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the single-clump candidates saw their masses increase slightly with the second approach. 
Although this is not a very significant effect, it may follow from the one element of the 
algorithm which changed, namely, that positions are fixed to the X-ray positions rather than 
marginalized over. The small difference between the two results implies that the X-ray peak 
position was a good proxy for the cluster center position, and that we did not marginalize 
over too large an area. 

The striking effect of clump multiplicity demonstrates that care must be taken not to 
introduce a systematic error in the measured cluster masses due to the presence of nearby 
shear peaks. It also implies that there could be significant systematic errors from nearby 
undetected masses. 

We now compare our results to previously published mass estimates. There is a dramatic 
change in Candidate 8. Recall the single-profile fit result of l. Tln? x lO^'^/i'^MfD pre sented 



in Table [2] was higher than, but consistent with, the corrected IWittman et al.l (120031 ) result 
of 2.8^1.0 X lO^^/i~^M0. When we simultaneously fit the additional extended X-ray source 
in the field, the additional source is assigned most of the mass, 2.4 ± 0.9 x 10^^/i^^Mq, and 
the main source is left with only 0.9 ± 0.5 x IO^^/i^^A/q. The reason is unclear: the main 
source has a larger X-ray flux by a factor of 1.4 despite being a t a higher redsh i ft (z = 0.68 



vs. 0.61). The dynamical and strong-lensing mass estimates in Wittman et al. (j2003 ) were 



also consistent with the higher mass, and cannot be affected by the secondary cluster at 
lower redshift. It could be argued that those estimates do not strongly rule out the lower 
mass estimate, given the small number of members with spectroscopy (17) and the unknown 
redshift of the strongly lensed arc. Hence, a complete understanding of this system may 
require more data, in the form of an arc redshift, Hubble Space Telescope imaging for a 
higher- resolution weak lensing analysis, and/or more spectroscopy. 

The sec o nd sy stem with previously published mass estimates is Abell 781, for which 
Sehgal et al.l (j2008l ) simultaneously fit NFW profiles to clumps la (which they called the 



"Main" cluster), lb (which they called the "Middle" cluster), and Ic (which they called the 
"East" cluster) o Our simultaneously-fit NFW masses are consistent with theirs. Although 
TableOgives the impression that the new results are systematically somewhat lower, the total 
discrepancy of the three clumps is only about one in terms of x^, for three degrees of freedom. 
This is not significant. The velocity dispersions inferred from the multi ple simultane o us SI S 



fits also agree with the spectroscopic velocity dispersions measured by iGeller et al.l (120051 ). 



^Note that ISehgal et al.l (j2008l ) considered cluster Id to be merging with cluster l a. and did not fit it 



separately. We do fit it separately, which should lower the mass of cluster la somewhat. Sehgal et al. ( 20081 ) 



also found a "West" cluster in XMM data covering a larger field. We do not consider that cluster here, as 
its low mass and large projected distance make it unlikely to affect the results. 
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When the shear fits are not simultaneous (Table [2]), the inferred velocity dispersions are 
significantly overstimated. 

Finally, we compare the results between the two different profiles, SIS and NFW, for all 
the measured masses (Figure H]). Only the masses of candidates la and lb from the multiple 
profile fit have more than a la discrepancy between the two profile types; their SIS masses are 
lower than their NFW masses by about 2cr. When the clumps are not fit simultaneously (blue 
points), masses inferred from the two profiles are in excellent agreement. Given the greater 
susceptibility of the SIS to boosting by neighbors, we expect a shift in this direction, but we 
would have expected the non-simultaneous fits to lie above the unity-slope relation rather 
than on it. Given that these are the two highest-shear clumps in the survey, this somewhat 
unexpected behavior may be indicating that high-mass neighbors induce additional subtle 
effects not accounted for by merely simultaneous fitting (for example, magnification changing 
the source redshift distribution). 

Figure m also suggests that the NFW profile tends to assign more mass to the candidate; 
nearly all the multiple-profile points lie below the line of unity slope. For the single-profile 
fits, this effect seems to be nearly cancelled by the tendency of the SIS to pull in more 
neighboring mass. We investigated this observation with data simulated using an NFW 
profile. We fitted the data with both SIS and NFW profiles and then compared the masses 
measured from both profiles. We found instead that fitting the SIS profile caused a slight 
over-estimation of the true mass, but this was insignificant given the size of the statistical 
errors. The over-estimation was due to the inner regions of the profile having more weight 
in the fitting, and in these regions the SIS density is less than the NFW density for the same 
mass. Therefore it is possible the reason why the NFW profile seems to assign more mass to 
the cluster than the SIS profile lies within the (unknown) exact details of the true density 
profile of clusters at these inner radii. This is reasonable because the shape of the cluster 
density profile at the lowest radii is the most uncertain. The NFW profile is predicted by 
cold dark matter only simulations, and would be affected at low radii by the baryonic cluster 
component, the existence of massive neutrinos or even the possibility that the dark matter 
has a non-zero cross-section. 



6. Summary and Discussion 

We have measured the masses of a sample of shear-selected clusters from the Deep Lens 
Survey by fitting NFW and SIS profiles to the full dataset of shears and source redshifts. 
For each type of profile, we fit the data in two different ways: (i) we fit to each shear peak 
separately, marginalizing over likely locations of the mass profile center; (ii) we fit to all the 
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extended X-ray sources simultaneously, fixing the mass profile locations at those of the X-ray 
peaks. We find that: 

• Masses measured with either profile are consistent with each other on a per-cluster 
basis. This is because the two profiles are similar to each other at the cluster radii 
probed by the weak lensing. There are hints that the NFW profile tends to assign more 
mass on average, but a larger and high er-S/N samp l e wou ld be required to confirm 



this. The same tendency was found by lOkabe et al.l (120091 ) . who studied 22 clusters 



with weak lensing signal measured from Subaru/Suprime-Cam imaging data. 

Fitting multiple clumps simultaneously can yield a much lower mass estimate compared 
to fitting them separately. This is because the single-profile fit confiates shear from 
nearby peaks with shear from the peak being fit. The size of this effect was a factor of 
about 1.5 to 2 for Abell 781 and for Candidate 2. 

Uncertainty in the mass profile center had a relatively small effect here, but should 
not be overlooked. For the candidates with a single extended X-ray source, the only 
difference between approaches (i) and (ii) was the marginalization over, vs. fixing of, 
the mass profile center, and the results did not change dramatically. There was a shift 
in the sense that marginalizing over possible centers did yield lower mass estimates, 
but it was small compared to the typical per-cluster uncertainty here. While adequate 
for the current purpose, with larger samples this may become a nontrivial systematic 
issue. Especially since we find from simplified simulations that a small offset of just 0.4 
arcminutes could cause a systematic bias on the mass if the precision of the statistical 
errors is greater than 5 per cent. More investigation into the typical BCG/X-ray and 
mass peak position offsets is required to make sure this systematic is fully under control. 

A few clusters have masses which are consistent or nearly consistent with zero. Al- 
though several effects may be contributing, the most likely is that the original shear 
selection, having been done without regard to source redshifts, was vulnerable to line- 
of-sight projections in a way that the current fitting tech nique is not. This pro blem 



may have been compounded by the heavy smoothing in the lWittman et al.l (I20061) con- 



yergen ce maps. As shown by the spectroscopy and X-ray detections in IWittman et al. 



(120061 ) ■ these are real clusters, but they may be so poor that they passed the selection 



only because of nearby projections. To mographic selection s hould reduce this prob- 



lem, but, according to the simulations of iHennawi fc Spergell (120051 ). cannot eliminate 



it. FoUowup spectroscopy and/or X-ray observations can prove that a cluster exists 
at that location, but much more extensive (perhaps prohibitive) foUowup would be 
required to prove that that are no additional groups projected nearby that boosted 
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the cluster into the sample. An alternative would be to think o f these syste n is sira ply 
as shear peaks rath er than clusters, and use the formalism of IWang et al.l (120091 ) or 



Marian et al.l (120081 ) to constrain cosmology. 



Comparing to previous mass estimates of some of these systems, we find that Abell 781 
is entirely consistent. However, Candidate 8 is more puzzling than previously thought; the 
new fits indicate that the secondary X-ray peak has more mass than the primary one. More 
data may be required to fully understand this system. 

Simulations suggest that the ellipticity of c lusters could be used as a cosrnologic al probe, 
if sufficiently good statistics can be obtained (IHopkins et al.ll2005l : IHo et al.ll2006l ). Obser- 
vations do suggest that cluster shapes, measured by the distribution of galaxies. X-rays, SZ 
effect and weak lensing, are usu ally not spherical. There is also indication that the mean 
ellipticity evolves with redshift (IMelott et al.l l200ll : iPlionid |2002| ). with clusters becoming 
less elliptical over time. This is expected because of hierarchical cluster formation and the 
reduced merger rate as clusters undergo a relaxation process. The systematic ellipticity of 
cluster dark matter h alos was recently detected by stacking hundreds of clusters in SDSS 
(lEvans fc Bridlelboooh . 



An extension to the current paper would be therefore, to allow the individual clumps 
to be fitted with an elliptical profile instead of an axi-symmetric one. This may reduce 
the number of clumps required in the fit, or produce more realistic mass estimates for the 
individual clumps. H owever the ellipticity of clusters is difficult to measure for individual 



20041 ) due to the high noise levels on individual shear 



systems (althou gh see ICvpriano et al. 

measurements. ICorless &: King] (120091 ) looked at the resulting biases on the cluster mass 
function from fitting a spherical model to both spherical and elliptical (triaxial) halos after 
stacking the lensing signals of many clusters binned by mass- correlated observables. They 
fitted a circular model to the stack of elliptical halos and measured M200 in a circular contour 
and found that this did not match with the M200 values of the original individual halos, if this 
M200 was calculated in elliptical isodensity contours. However they found that the stacked 
M200 matched well to the M200 of the individual halos, if the M200 values of the original 
elliptical halos were defined using circular contours instead, where the density at different 
points around these circular contours was not constant. 

We performed a simulation to test this effect on a per-cluster basis. The clusters in 
Corless fc King] (120091 ) were stacked without attention to aligning their axes, so the random 
orientations of many halos in each bin would artificially circularise the average of the stacked 
halos. In our simulations we calculate the difference between (i) the mass obtained from 
fitting a circular model and (ii) the mass within the circle containing 200 pcrit of the best-fit 
elliptical model. We find these two numbers are the same to better than 1 per cent accuracy. 
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Therefore we conclude that there is neghgible systematic bias in our results from fitting 
a spherical profile to an elliptical halo as long as M200 in circular radii is the quantity of 
interest . 

Many different cluster surveys will be taking place in the near future. They will aim 
to probe the growth of structure over time, beginning a new genre of cosmological analysis. 
Understanding the sample selection and relating the cluster observables to the mass is a 
large task. Shear selection in principle eases this task because it is independent of the nature 
of the matter and of its dynamical state, but there are clearly several issues which must be 
resolved, principally the observationally difficult task of accounting for projections, but also 
including determination of the mass profile center without biasing the mass estimate. In 
the future, large surveys such as the Large-aperture Synoptic Survey T elescope (LSST) will 
find on the order of 100,000 shear-selected clusters (iTyson et al.ll2003l ). With such massive 
surveys systemat ics will become important. The modest samples produced by DLS a nd other 
current surveys (IMiyazaki et al.l 120071 : iGavazzi fc Soucaill 120071 : [Dietrich et al.l 120071 ) will be 
the testing ground for developing methods to successfully analyze these new large data sets. 
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